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ABSTRACT 

We present a luminosity function for low surface brightness (LSB) galaxies 
identified in the APM survey of [Impey et ah (1996)| . These galaxies have 
central surface brightnesses (p(0)) in B in the range 22.0 < mu{0) < 25.0. 
Using standard maximum-likelihood estimators, we determine that the best-fit 
Schechter function parameters for this luminosity function (LF) are a = —1.42, 
M* = —18.34, and (j)* = 0.0036, assuming Rq = 100 hioo kms“^ Mpc“^. We 
compare the luminosity and number densities derived from this luminosity 
function to those obtained from other recent held galaxy studies and hnd that 
surveys which do not take account of the observation selection bias imposed by 
surface brightness are missing a substantial fraction of the galaxies in the local 
universe. Under our most conservative estimates, our derivation of the LF for 
LSB galaxies suggests that the GfA redshift survey has missed at least one third 
of the local galaxy population. This overlooked fraction is not enough by itself 
to explain the large number of faint blue galaxies observed at moderate redshift 
under no-evolution models, but it does help close the gap between local and 
moderate-redshift galaxy counts. 
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1 Introduction 


The optical luminosity function (LF) of galaxies is one of the fundamental building blocks 
of cosmology. Accurate knowledge of the luminosity function is necessary for, among other 
things, estimating the mean luminosity density of the universe, and predicting the redshift 
distribution of objects in various magnitude intervals (see e.g., the review by Binggeli et ah 
1988). The shape of the luminosity function also provides an important test for theories of 
galaxy formation (e.g.. Press & Schechter 1974). Further, considerable attention has been 
focussed of late on the large numbers of blue galaxies found in deep surveys, hrst described by 
Kron (1980)1 and [Hall fc Mackay (1984)| . The degree to which number counts of these galaxies 
exceed those predicted from local observations (e.g., Bruzual & Kron 1980 and Guiderdoni 
&: Rocca-Volmerange 1990), and indeed whether an excess exists at all (compare Koo et al. 
1993 and McGaugh 1994), depend on the shape, normalization and color dependence of the 
luminosity function. 

One of the problems with building a galaxy luminosity function is that surveys are 
limited in the detection of diffuse galaxies by the brightness of the night sky, and in the 
detection of compact galaxies by the difficulty in distinguishing stars and galaxies. As 
Disney (1976) and Disney &: Phillipps (1983] have demonstrated, at a given luminosity a 
survey will identify preferentially those galaxies that have the maximum possible angular size 
above the limiting isophote. At a constant luminosity, galaxies of high surface brightness 
(HSB) become indistinguishable from stars, and galaxies of low surface brightness (LSB) 
fall below the limiting isophote over most of their extent. Although they purport to be 
magnitude limited, galaxy surveys which do not take account of surface brightness effects 
are missing an unknown but potentially large number of galaxies in each magnitude bin. 
Recent surveys of the Virgo cluster by [Impey et al. (198K]| and of the Fornax cluster by 
Irwin et al. (199'0]| and [Bothun et al. (1991)| have taken account of this potential source 
of bias by deliberately searching for LSB galaxies. They have found that previous surveys 
missed a signihcant fraction of the cluster populations, particularly at fainter luminosities 
{Mb ^ —16), and [Impey et al. (198'^ determined that inclusion of LSB galaxies in Virgo 
steepened the low-luminosity tail of that cluster’s luminosity function considerably. To date, 
however, no estimates of the held galaxy luminosity function have addressed the effects of 
surface brightness bias. However, [McGaugh et al. (1995a)| found that the space density of 
galaxies as a function of central surface brightness appears to be hat below /iB(0) = 22.0. 
Also, [Sprayberry et al. (1996)[ found a space density of galaxies as a function of central 
surface brightness that appeared hat below //^(O) = 23.0 after descending from a peak 
around /iB(0) = 21.75. Although many of these LSB galaxies are not necessarily faint, the 
forms of these distribution functions strongly suggest that the normalization of the galaxy 
space density at 2 ; = 0 has been strongly inhuenced by surface brightness selection ehects. 


We have recently completed a survey for LSB galaxies in the region dehned by 
—3° < S <3° and |6| > 30°, surveying about 786 square degrees of sky with the Automated 
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Plate Measuring (APM) system at Cambridge.^ We have identified 693 galaxies, most 
previously uncataloged and most with central surface brightness > 22 magarcsec”^. 

The complete catalog of this survey appears in |lmpey et ah (1996]| (Paper I). The selection 
effects and completeness corrections for the survey are analyzed in detail in ^prayberry et al. 
(1996)1 (Paper II). 


In this paper, we present the luminosity function for LSB galaxies from the APM survey 
and compare that luminosity function to those obtained from the CfA redshift survey. We 
also review suggestions by [Phillipps et al. (1990)| , [McGaugh (1994)| , [McLeod (1994)| , and 
Ferguson fc McGaugh (199F)| that LSB galaxies might account at least partially for the large 
numbers of faint blue galaxies seen in deep surveys. Section describes the survey data and 
presents the samples used for determining the luminosity function and the corrections applied 
to those samples. Section ^ covers the methods used to develop the luminosity functions. 
Section ^ presents the luminosity functions and compares the results to those obtained from 
the CfA redshift survey. Section reviews the consequences of this LSB luminosity function 
for the general field luminosity function and for the question of local counterparts to the 
faint blue galaxies. Finally, Section summarizes our conclusions. Throughout this paper, 
we assume Hq = 100 hiookms"^ Mpc“h Also, all magnitudes and surface brightnesses used 
here are in the Johnson B band. 


2 Samples Used 


The APM survey for LSB galaxies is presented in Paper I, and Paper II describes the details 
of how LSB galaxies were identified and calibrated. Paper II also presents a selection function 
that gives the completeness of the survey as a function of galaxy central surface brightness 
and scale length (hereafter, “the APM selection function”). 


We conducted followup optical spectroscopy at the Multiple Mirror Telescope^ and 21 
cm H I spectroscopy at Arecibo Observatory^ to obtain radial velocities for as many of 
the galaxies as possible. To date we have measured recessional velocities for 332 of the 
693 galaxies on the list, of which 190 come from H I spectroscopy and 142 from optical 
spectroscopy. These heliocentric velocities are presented in Paper 1. For developing the 
luminosity function, we have further corrected these heliocentric velocities to the rest frame 
of the Local Group, using the standard correction Vcorr = "yhei + 300sin/ cos6. No correction 
was applied for Virgocentric infall since the median velocity of the sample places most of 


^The APM is a National Astronomy Facility, at the Institute of Astronomy, operated by the Royal 


Greenwich Observatory. A general description of the APM facility is given by Kibblewhite et al. (1984). 

^The Multiple Mirror Telescope is a facility jointly operated by the Smithsonian Institution and the 
University of Arizona. 

^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center. The NAIC is 
operated by Cornell University under a cooperative agreement with the National Science Foundation. 
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Figure 1: Structural properties of the complete LSB sample and the subset with radial 
velocities, (a) shows the distribution as a function of B central surface brightness, and 
(b) shows the distribution as a function of half-light radius. In the upper panels, the 
dotted histogram is the distribution of the complete sample, and the solid histogram is 
the distribution of the subset with velocities. In the lower panels, the filled circles show the 
fraction of galaxies with velocities for each bin, with error bars from counting statistics. The 
solid lines show the parametrizations described in the text. 

the galaxies well beyond the Local Supercluster. These corrected velocities were then used 
to estimate distance moduli using the relation: 

m- M = 5 [log Vcorr “ log i^o + 5] (1) 

assuming as noted above that Hq = 100 hiookms“^ Mpc“^. 

The galaxies with velocities do not form a random subset of the overall survey. For 
reasons of observational efficiency, like all other galaxy surveyers we favored galaxies of 
higher central surface brightness and larger angular size. Figure shows the distributions of 
central surface brightness and half-light radius for the complete sample and for the subset 
with velocities, along with the ratios of the two sets by bin. We assume that the galaxies 
for which we have measured redshifts are representative of all galaxies in a given bin of 
surface brightness and angular size. This additional source of bias must be taken into 
account in preparing a luminosity function. We have parameterized this bias in the simple 
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forms depicted in Figure ^ three separate linear fits in the different regions of the /i(0) 
distribution 

r 1.000, /i(0) < 20.25 

i 4.950 -0.194/i(0), 20.25 < p(0) < 25.0 (2) 

[ 0.111, /i(0) > 25.0 

where /i(0) is in magarcsec”^, and in the different regions of the half-light radius distribution 

r 0.667, Teff < 3 

= S —0.130 -|- 0.076 Te//, 3 < Te// < 13 (3) 

[ 0.773, r,ff > 13 

where rg// is in arcseconds. The hnal probability that an LSB galaxy will be detected by 
the APM and included in the subset with velocities is given by 


Ptot = PAPM X P^ X Pr^ 


( 4 ) 


where Papm is the probability derived from the APM selection function of Paper II. 
Equation ^ assumes that the corrections in p(0) and r^ff are separable. This assumption is 
reasonable for our sample, because p(0) and r^ff are uncorrelated: Pearson’s r = 0.075 and 
the Spearman rank correlation coefficient s = —0.121, and neither coefficient is signihcantly 
different from zero. 


We note that Figure |l| shows that the surface brightness range 23.5 < p(0) < 24.5 
includes a large number of identihed galaxies, but that a very small fraction of those galaxies 
were observed spectroscopically. Also, the observed fraction as a function of angular size 
declines sharply at small sizes. These features are artifacts of the two stages in which 
the APM survey was performed. The hrst stage identihed LSB galaxies of large angular 
size, and all the followup spectroscopy was performed on galaxies in this hrst list. The 
second stage identihed small angular size galaxies, which also tended to be predominantly 
in the surface brightness range 23.5 < ju(0) < 24.5. The interested reader is referred to 
Paper II for a more complete discussion of the survey mechanics. Here we note only that the 
actual observed fraction in the range 23.5 < ju(0) < 24.5 lies below the parametrization of 
Equation which implies that the parametrized correction is too small for those two surface 
brightness bins. Any bias introduced by this ehect is “conservative”, in that it will result in 
an underestimation of the total number of LSB galaxies. 


We can estimate the completeness of our sample of galaxies using the (V/Vmax) test of 
^chmidt (1968)| . For the complete set of 693 galaxies identihed by the APM, the test yields 
(V/Vmax) = 0.15 ± 0.04 with no corrections for incompleteness, and (V/Vmax) = 0.44 ± 0.06 
after correcting for incompleteness using the APM selection function described in Paper IF 
For the subset of 332 galaxies with velocities, the test gives (V/Vmax) = 0.04 ± 0.05 with no 
corrections for incompleteness, (V/Vmax) = 0.34±0.07 after applying just the APM selection 
function, and (V/Vmax) = 0.50 ± 0.07 after applying the APM selection function and the 
further correction for incompleteness in the velocity observations from Equations 0, 0, and 0 
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(as depicted in Figure Q). The corrections thus substantially remove the incompleteness in 
both the complete set and in the subset chosen for spectroscopy. 


There is yet another source of bias to be found in the magnitudes measured for LSB 
galaxies. The magnitudes measured in our survey are isophotal magnitudes, not extrapolated 
or asymptotic. The median limiting isophote is fium ~ 27.4 magarcsec”^. As authors from 
pisney (1976)| to [McGaugh (1994)| have pointed out, use of isophotal magnitudes will cause 
galaxy luminosities to be underestimated, and the underestimation becomes more severe with 
decreasing central surface brightness. Most LSB galaxies are well-described by exponential 
surface brightness prohles (Impey et ah 1988, Bothun et ah 1991, and McGaugh & Bothun 
1994) of the form 

/i(r) = /i(0) + 1.086y (5) 

where p(0) is the central surface brightness in magarcsec“^ and I is the exponential scale 
length in arcseconds. This simple analytical form allows a direct calculation of the ratio of 
the total galaxy flux to that observed within the limiting isophote, as 


F, 


obs 


F 


= 1 - {1 + ni)e-^^ 


( 6 ) 


tot 


where ni is the number of scale lengths I observed within the limiting isophote. This simple 
approximation will clearly understate the ratio for galaxies with central condensations, such 
as spirals with bulges. The isophotal aperture in units of the galaxy scale length is then 
given by 

h/im - h(0) - 10 log(l + z) - k{z) 


ni = 


1.086 


( 7 ) 


where is the surface brightness of the limiting isophote. The hrst term involving z 
accounts for the (l-f-z)"^ cosmological dimming in surface brightness, and the second corrects 
for the redshifting of the galaxy’s spectral energy distribution (the k correction). The k 
correction of course depends on galaxy type as well as redshift. The magnitudes and surface 
brightness for the LSB galaxies with velocities have been corrected as described in Paper II 
using the tabulated k corrections of Coleman et ah (1980) . The B — V and V — R colors for 
galaxy types Sbc, Scd, and Irr closely match the range of colors observed among the galaxies 
for which we obtained CCD photometry. The absolute magnitudes have been corrected 
according to Equations and so as to avoid skewing the luminosity function by this 
tendency to underestimate galaxy luminosities. 


Of course, our set of LSB galaxies is not itself a fair sample of the local galaxy population, 
precisely because it excludes most galaxies with /r(0) ^ 22 magarcsec”^. However, it is still 
useful to derive a luminosity function for this set, so that this LF can be compared to one 
derived from higher surface brightness galaxies. In this way, it is possible to obtain some 
idea of how surface brightness selection effects have influenced estimates of the density of 
local galaxies (see also McGaugh et ah 1995 and Paper II). To validate such a comparison, 
it is necessary hrst to compare the range of surface brightnesses covered by the present set 
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of LSB galaxies with the range covered by other surveys. Unfortunately, no other recent 
galaxy redshift surveys have published surface brightness data for their galaxies. Thanks 
to the recent release of a digitized version of the original Palomar Observatory Sky Survey 
(the Digitized Sky Survey^ or DSS), it is now possible to make independent measurements 
of the basic photometric parameters of any object visible on the original survey, when the 
celestial coordinates of the surveyed galaxies are known. The CfA Redshift Survey described 
by e.g., [IVfarzke et al. (1994b)| is based on Zwicky’s Catalog of Galaxies and Clusters of 
Galaxies, which was in turn created by visual examination of the Palomar Observatory Sky 
Survey plates, so every object included in that survey should be visible on the DSS. Most 
importantly, the coordinates of galaxies surveyed by the CfA are publicly available, so that 
it is possible to retrieve images of the surveyed galaxies from the DSS. Thus it should be 
possible to measure the surface brightness range covered by the CfA Redshift survey. The 
lack of publicly available coordinates prevents us from making a similar analysis of other 
recent redshift surveys. 

We recovered from the Astrophysics Data System listing of the CfA Redshift Survey 
the coordinates of every galaxy listed in the regions of sky used by [Marzke et al. (1994b )| . 
We subdivided that list according to the morphological categories used by [Marzke et al. 
|(1994a)| , and we randomly selected 10% of the galaxies within each morphological class to 
keep the number of galaxies manageable. This selection yielded a list of 579 galaxies. We 
then retrieved images from the DSS of this randomly chosen subset and analyzed the images 
using the same algorithms used in our APM LSB galaxy survey. In this way, we obtained 
extrapolated central surface brightnesses for the CfA galaxies that are directly comparable 
to those obtained in the course of the APM survey. Paper II contains a complete description 
of the process of estimating the extrapolated central surface brightness. As a check on the 
calibrations, we also retrieved from the DSS images of a randomly chosen subset of the 
APM LSB galaxies and analyzed them. After cross-calibration, the results for the APM 
LSB galaxies were consistent with those obtained from the deeper UKST plate materials 
used in the APM LSB survey, with the exception that the lowest surface brightness objects 
were not visible on the DSS. 


The surface brightness distribution for the CfA Redshift survey is shown in the upper 
panel of Figure ^ The solidly drawn smooth curve represents the best Gaussian fit to the 
GfA distribution. The lower panel shows the complete SB distribution obtained by the 
APM for one UKST held. Also drawn for illustration in each panel is a dashed curve 
representing the canonical “Freeman Law”, a Gaussian centered at /i(0) = 21.65 with 
a = 0.35 (Freeman 1970). It is clear from Figure that the range of surface brightnesses 
covered by the GfA Redshift Survey is very narrow, narrower even than the “Freeman Law.” 

^Based on photographic data of the National Geographic Society - Palomar Observatory Sky Survey 
(NGS-POSS) obtained using the Oschin Telescope on Palomar Mountain. The NGS-POSS was funded by 
a grant from the National Geographic Society to the Galifornia Institute of Technology. The plates were 
processed into the present compressed digital form with their permission. The Digitized Sky Survey was 
produced at the Space Telescope Science Institute under US Government grant NAG W-2I66. 
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Figure 2: Distributions of central surface brightness for (a) a randomly chosen sample of 
galaxies from the CfA redshift survey, and (b) the complete list of galaxies identihed by 
machine scan of one UKST survey held. In (a) the solid curve represents the best ht of 
a Gaussian to the CfA survey surface brightness distribution. In both panels the dashed 
Gaussians illustrate the canonical “Freeman Law” of /i(0) = 21.65 ± 0.35. The distribution 
in (b) is corrected for incompleteness of the detection algorithm for /^^(O) ^ 25 as described 
in Paper II. 


The best-ht Gaussian to the CfA distribution has a center at /i(0) = 21.44 and a = 0.19. 
This is completely consistent with the investigation of the Zwicky magnitude scale by [Bothun 
&; Cornell (1990)| who hnd that this magnitude is not a sky-limited magnitude. In this case, 
one expects surface brightness effects to completely dominate the magnitude estimates. In 
essence, the Zwicky magnitude is very much a “bulge” or high surface brightness magnitude 
and is insensitive to extended, low surface brightness light. In contrast, the APM LSB 
survey has identihed galaxies over a much broader range, as described in Paper II. Clearly, 
the identihcation of galaxies for the CfA Redshift survey suffered from a substantial bias 
against LSB galaxies. In all the following analysis, we use only those galaxies from the APM 
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M, 


Figure 3: Distribution of absolute magnitudes for the LSB galaxies (/i(0) > 22.0 
magarcsec“^) used to develop the LF. This distribution includes the effects of the correction 
from isophotal to total magnitudes described in Equations and 


LSB survey with /i(0) > 22.0 magarcsec“^, or 3 ct fainter than the typical value found in 
the CfA Survey. This limitation assures that the resulting LF covers a different regime of 
surface brightness parameter space from that covered by the LFs of |Marzke et ah (1994b) 
and [Marzke et ah (1994a)| . We note that there is a weak LSB tail in the CfA distribution: 
the overall xt Gaussian £t is 1.37, virtually all of which is due to this tail. However, 

the very weakness of this tail, when compared to the APM distribution in the lower panel, 
underscores the severity of the SB selection bias inherent in the CfA survey. We note also 
that the CfA survey does not identify nearly as many high surface brightness galaxies as 
does the APM. This lack is most likely due to the general absence of galaxies smaller than 1 
arcminute from the Zwicky catalog; many of the high surface brightness galaxies identified 
by the APM are smaller than 1 arcminute. Figure ^ shows the distribution of absolute 
magnitudes for the LSB survey galaxies with /i(0) > 22.0 magarcsec”^. 


3 Methods 

The differential luminosity function of held galaxies (j){M)dM is defined as the function 
giving, at each absolute magnitude M, the number of galaxies per Mpc“^ in the luminosity 
interval M + dM/2 < M < M — dM/2. Because the area surveyed by the APM LSB 
survey covers a wide area of sky and cuts across several large scale structures, we adopt two 
density-independent techniques for estimating the LF. The hrst is the parametric maximum 
likelihood technique developed by jSandage et ah (1979)| (hereafter STY). The second is 
the stepwise maximum likelihood method (hereafter SWML) developed by [Efstathiou et al. 
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(1988). Both methods assume that the LF has a universal form, independent of position, 
allowing the probability of a galaxy’s inclusion in a complete catalog to be written simply in 
terms of the LF itself. The STY method is continuous and uses all the galaxy data, but it 
requires the assumption of a parametrized form for the LF. It therefore gives no information 
as to the suitability of the parametrized form chosen to represent the LF. The SWML 
method requires binning the data, but it requires no assumptions about the shape of the LF. 
It can therefore be used in combination with the STY method to provide an independent 
check on the goodness-of-£t of the chosen parametrization, as described by [Etstathiou et al.| 
(1988 )|. Like [IVfarzke et al. (1994a)| and virtually all others who have used this combination 


of methods, we assume in the STY method a luminosity function parameterization in the 
form hrst proposed by ^chechter (1976) , which is written in absolute magnitudes as 


0(M)dM = O.41nlO0* (10 




-M)) 


dM 


( 8 ) 


Using the two methods together thus gives best-£t values for the Schechter function 
parameters a (the faint-end slope) and M* (the characteristic absolute magnitude of the 
“knee”), as well as a probability that the underlying galaxy population is well-described by 
the best-£t Schechter function. 

There is one major difhculty with applying these methods to the APM LSB galaxy 
survey data. Both the STY method and the SWML method assume that the galaxy catalog 
in use is magnitude limited, or that all galaxies with m < mum have the same probability 
{p < 1) of being included in the catalog, as in the case of a redshift survey that uniformly 
samples a magnitude-limited catalog with 1/n sampling. In our case, however, each galaxy 
has a unique probability of inclusion that is determined from Equation so the given forms 
of the STY and SWML methods require modihcation. ^ucca et al. (1994]| recently addressed 
this problem. They derived a simple modihcation to the STY estimator that accounts for 
the unique observation probability assigned to each galaxy: 


N 

^ = t[p7 

i=l 


(9) 


where C is the likelihood to be maximized, the weight Wi is defined as the inverse of 
the probability that the ith galaxy will be included in the sample (i.e., for our situation 
Wi = l/ptot,h with ptot,i from Equation H), and pi is as dehned by STY: 


Pi = (f){Mi) 



(i){M)dM 


( 10 ) 


maxi^Zj^ 


The corresponding change to the SWML estimator of [Efstathiou et al. (1988)| immediately 
yields: 

N N ( Np 'I 

In £ = ^ lU {Mi - Mk)wi In ^ Wi In < X] {Mmax{zi) - Mj^ > -F const (11) 


2 = 1 


2 = 1 


J = '^ 
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where the (pk are the luminosity function values within each bin, N is the total number 
of galaxies in the sample, Np is the number of steps, Mmax{zi) is the maximum (i.e., the 
faintest) absolute magnitude visible at Zi, AM is the bin width in magnitudes, and the 
window functions are 


W{x) 


1, |a;| < AM/2 

0, otherwise 


( 12 ) 


and 

r 0, a; < -AM/2 

H{x) = l {x/AM+ 1/2), \x\<AM/2 

[ 1, x> AM/2 

There is an implied sum over the doubled index k in the hrst term of Equation m 


(13) 


Finally, the survey biases must also be incorporated into the normalization. Both the 
STY and SWML estimators are normalized in the manner described by |Etstathiou et al. 


1(1988)1 using the unbiased minimum variance estimate of the mean density as developed 
by 0 avis fc Huchra (1982)| , but with a modihcation to the estimator to incorporate the 
corrections for survey incompleteness. This normalization proceeds in three steps. First, a 
selection function is dehned as 


pMi j pMi 

S{x) = / (j){M)dM / / (j){M)dM (14) 

/ J M2 


for galaxies in the range Mi < M < M 2 , where M^axix) is the maximum (i.e., the faintest) 
absolute magnitude visible at distance x according to the catalog limits. Second, this 
selection function is then corrected to incorporate the incompleteness correction, so that 
it includes the combined probability of detecting and spectroscopically observing an LSB 
galaxy in our survey: 

Stoti+i) S{x/) X Ptot,i (1^) 

where Ptot,i is obtained from |^. Finally, the mean density of galaxies is obtained from the 
corrected selection function as described by [Etstathiou et al. (1988)| : 


(n) 


N 


-E 



(16) 


where the sum extends over all the galaxies in volume V. The mean density is converted to 
a Schechter function normalization as: 


_M_ 

T (a + 1, 100'4 (m*-M2)) _Y(a + l, 


(17) 


where T is the Euler incomplete gamma function. 

I ' J 

Zucca et al. (1994) also estimated the effects of failing to consider the individual galaxy 
weights. Their simulations revealed that use of Equation |T^ to determine the Schechter 
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Figure 4: Luminosity function for LSB galaxies from the APM survey. The solid line 
represents the maximum likelihood Schechter function, and the points with error bars 
represent the model-independent step-wise maximum likelihood function. Note that for the 
LSB galaxies, the model-independent binned LF shows a signihcant excess of low luminosity 
galaxies, beyond the level of the maximum likelihood Schechter function. The dashed line 
shows the maximum likelihood Schechter function estimated by pVlarzke et al. (1994a)| for 
all morphological types in the CfA redshift survey, and the dotted line shows the maximum 
likelihood Schechter function estimated by Marzke et al. for irregulars in the CfA redshift 
survey 


function parameters for a galaxy sample with signihcant incompleteness {{V/Vmax) ^ 0.3) 
would bias the results towards hatter faint-end slopes (i.e., lower absolute values of a) 
and brighter values of M*. We can objectively determine individual galaxy weights from 
parameters of our survey technique (the APM selection function) and from the internal 
statistics of our followup observations (Figure |I] and Equations || and so Equations ^ and 
0 are the clear techniques of choice for our data. 


4 Results 

Figure shows the luminosity function for the LSB galaxies (/i(0) > 22.0) from the APM 
survey. The solid line represents the maximum likelihood Schechter function from the STY 
method, and the points with error bars represent the model-independent SWML method. 
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Figure 5: Histogram of T-types for LSB galaxies from the APM survey. Galaxies that 
appeared to be spirals but whose images on the APM scans were too small to permit reliable 
further classification were assigned T = 5, so the number in that bin is somewhat inflated. 
Interacting galaxies were assigned T = 11. 


As is obvious from Figure ^ the maximum likelihood Schechter function is a very poor 
representation of the “true” distribution as determined by the SWML method: the reduced 
xl from the likelihood ratio test of [Efstathiou et al. (1988)| is 14.04, which implies that 
the probability of exceeding this xl ^>7 chance is ~ 2.5 x 10“^^. The Schechter function is 
particularly poor at the low-luminosity end. There, the model-independent SWML method 
hnds two to three times the galaxy density predicted by the maximum likelihood Schechter 
function. The SWML bins at M = —16, —15, —14, and —13 contain 31, 9, 15, and 12 
galaxies respectively. Across those four bins, the median correction due to the APM selection 
function (Paper II) is 0.759, and the median correction due to the incomplete spectroscopic 
observations (Equations |] and |^) is 0.259; the median total incompleteness correction 
(Equation is therefore 0.197. Our sampling of galaxies in these low-luminosity bins is 
quite sparse, and hence the uncertainties at this end are large. The formal uncertainties 
shown in Figure ^ may well understate the true range. 


The dashed line in Figure ^ represents the Schechter function estimated by |Marzke et al. 
(1994a)| for all galaxy morphologies in the CfA Redshift Survey. At the faintest luminosities, 
LSB galaxies in the range 22.0 < p(0) < 25.0 are more numerous than the HSB galaxies 
sampled by the CfA Survey if the comparison is based on the model-independent SWML 
points for the LSB galaxies, or approximately as numerous if the comparison is based on the 
maximum likelihood Schechter function. For all galaxies brighter than Mb < —15 the HSB 
galaxies are significantly more numerous. Table |l] lists the maximum likelihood Schechter 
function parameters for the LSB galaxies in the present study, along with similar parameters 
for all morphological types and for irregular galaxies from the CfA Survey. 
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The LSB sample from the APM survey is not restricted as to morphological type. It 
includes a few dwarf ellipticals and early spiral types. However, it is dominated by very late- 
type spirals and irregulars. As Figure |] shows, over half the LSB sample have de Vaucouleurs 
T-types of 9 or 10. For that reason, we have also shown the Schechter function derived by 

for irregulars (which they dehne as 8 < T < 10) in Figure ^ and 
Table |l|. The Schechter function for the CfA irregulars bears a striking resemblance to that 
derived here for the LSB galaxies. The steep low luminosity tail of the function for CfA 
irregulars seems to match the model independent SWML points for the LSB galaxies quite 
well. This similarity in LF slopes could be used as an argument that the high space density 
of star forming irregulars are the parent population of the fainter LSBs. Unfortunately, 
photometric surveys of LSBs continue to hnd no relation between SB and color which is 
required to support such a fading model. 


Marzke et ah (1994a) 


Table 1: Comparison of Luminosity Function Model Parameters 


Model/Survey 

a 

M* 

0* 

(1) 

(2) 

(3) 

(4) 

Maximum-Likelihood Schechter Functions: 
LSB 

-1.46 

-18.66 

0.0036 

CfA (all types) 

-1.02 

-18.90 

0.0201 

CfA (Sm - Im) 

-1.87 

-18.79 

0.0006 

Schechter Function -|- Power Law: 

LSB (giants) 

-0.92 

-18.19 

0.0060 

LSB (dwarfs) 

-2.20 

-16.00 

0.0041 


Notes: M* in B mag, 0* in /i^qq Mpc ^mag L 


Several authors have suggested that the LF for faint galaxies may exhibit an upturn from 
the pure Schechter form at faint luminosities. The LF of |lmpey et ah (1988)| for LSB dE’s in 
the Virgo cluster turns up at an apparent magnitude = 17; the increase is so steep that 
they were unable to rule out a divergent faint-end slope (i.e., a = —2.0). Upturns from the 
Schechter form have been observed in Coma by [Thompson fc Gregory (1993)| ; in nearby local 
groups by [Ferguson fc Sandage (1991)]; in four local Abell clusters by pe Propris et al. (1995) 


and in Coma, Abell 2554 and Abell 963 by priver fc Phillipps (1996) 


These deviations 
n < Mb < -15, 
In contrast. 


from the Schechter form are generally seen to begin in the range 
after adjustment to the distance scale used here {Hq = 100 hiookms"^ Mpc“^ 

[Ferguson fc Sandage (1988)[ found an LF for the Fornax cluster that was consistent with a 
single Schechter function having a faint-end slope of a = —1.34. We note that the form of 
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M, 


Figure 6: Luminosity function for LSB galaxies from the APM survey, determined using 
the two-component model described in the text (dashed curve), and maximum likelihood 
Schechter function (solid curve). The points with error bars represent the model-independent 
step-wise maximum likelihood LF. For galaxies brighter than M = —16 the two-component 
model is a standard Schechter function with all three parameters allowed to vary. For galaxies 
fainter than M = —16, the model is a power law with the normalization at M = —16 
constrained to match the value of the Schechter component there. Note that this model was 
£t to the SWML binned LF, and not to the individual galaxy data points. The hducial 
line labeled “Divergent” illustrates the faint end slope a = —2 where the integral of the LF 
becomes divergent. 


the SWML data points in Figure ^ is consistent with earlier hndings of a sharp change in 
slope at faint luminosities: the binned model-independent data points clearly break up from 
a smooth Schechter form at M = —16. 

To investigate this break in more detail, we also fit a two-component model to the SWML 
LF representation. For absolute magnitudes M < —16 (“giants”), the model followed the 
usual Schechter function form, with all three parameters (a, M*, and 0*) allowed to vary. 
For absolute magnitudes M > —16 (“dwarfs”), the model followed a simple power law, 
with the normalization at M = —16 constrained to match the Schechter model value there 
(i.e., only the slope a was allowed to vary in the £t). Results of this £t are depicted in 
Figure P, and the htted model parameters are listed in Table |l|. Two aspects of this £t 
deserve special comment. First, it is not possible to compare directly the goodness-of-£t of 
this two-component model to that of the STY maximum-likelihood Schechter function. The 
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Figure 7: Cumulative luminosity functions for LSB galaxies and for galaxies from the 
CfA redshift survey, plotted linearly as a function of Log L. The top axis represents the 
corresponding B magnitude. The cumulation runs from high to low luminosities (i.e., from 
right to left). The points represent the binned model independent LF for the LSB galaxies, 
the solid line represents the maximum likelihood Schechter function for the LSB galaxies, 
and the long-dashed line represents the two-component model LF for the LSB galaxies. 
The short-dashed line is the Schechter function determined by Marzke et ah (1994) for all 
morphological types in the CfA redshift survey, and the dotted line is the Schechter function 
determined by Marzke et al. (1994) for irregulars. 


quoted above for the STY result is obtained by a likelihood-ratio test which, like the 
STY model itself, is computed from the individual galaxy data points. The two-component 
model is a £t to the binned SWML LF, not to the individual galaxy data points, and thus 
its much larger xl {xt = 294.8, due to the much smaller number of degrees of freedom) is 
computed in a very different manner. Second, the two-component model yields results which 
are not physically plausible. The resulting slope of the dwarf galaxy power law is a = —2.20, 
which implies an infinite total luminosity if the two-component LF is integrated from zero 
to inhnity. This faint-end slope is best interpreted as a finding that the model-independent 
SWML indicates a very steeply increasing density of faint galaxies, so steep in fact that a 
divergent LF cannot be ruled out (cf. Impey, Bothun & Malin 1988). 


Because the low luminosity tail of the LSB luminosity function rises so steeply, the 
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contribution of LSB galaxies to the overall number density of field galaxies locally is 
quite large. Figure ^ shows the same luminosity functions as Figure ^ along with the 
two-component LF described above, but cumulated to show total number densities. The 
cumulation runs from high to low luminosities (i.e., from right to left), and the vertical axis 
scaling is linear. Across the whole range of luminosities, the LSB galaxies are almost twice as 
numerous as all the HSB galaxies in the CfA survey if the comparison is based on the SWML 
points for the LSB galaxies, or half as numerous if the comparison is based on the maximum 
likelihood Schechter function. Thus, even under the most conservative estimate, surveys like 
the CfA redshift survey have missed at least one-third of the local galaxy population due 
to surface brightness selection biases. The true missed fraction is almost certainly higher, 
even by the most conservative estimator. The LSB LF presented here covers only the range 
22.0 < /r(0) ^ 25.0, but as [McGaugh et al. (1995a)| and Paper II showed, the distribution 
appears flat for SB levels /i(0) > 25.0. Thus, surveys sensitive to fainter SB levels should 
hnd even higher number densities of galaxies. 


Despite the signihcance of their total numbers, LSB galaxies contribute little to the 
total luminosity density of the local universe, because the highest number densities of LSB 
galaxies occur at the lowest galaxy luminosities. Figure ^ shows the cumulative luminosity 
densities for LSB galaxies from the APM survey and HSB galaxies from the CfA redshift 
survey. Again, the cumulation runs from high to low luminosities (right to left) and the 
vertical scaling is linear. By any estimator, LSB galaxies represent a small fraction of the 
total luminosity emitted by HSB galaxies from the CfA survey. 


5 Implications 

Much attention has been devoted over the past 15 years to the population of faint blue 
galaxies revealed in surveys sensitive to extended objects as faint as uiej ~ 27. These 
galaxies become bluer at fainter apparent magnitudes (Lilly et al. 1995). They generally are 
not at extreme redshifts: [Lilly et al. (1995]| found a median redshift of Zmed ~ 0.56 for a 
sample of galaxies in the magnitude range 17.5 < Iab < 22.5. These galaxies are clustered 
more weakly than are most local bright galaxies, though their clustering strength is roughly 
comparable to that of local galaxies undergoing rapid star formation, per [Bernstein et al.j 
[( 1994)[ . Their numbers are signihcantly in excess of expectations based on local galaxy 
populations in the absence of evolution (Tyson 1988; Lilly et al. 1991; McLeod & Rieke 
1995). This excess has led some authors to suggest non-standard cosmologies as a possible 
explanation (Yoshii 1993), and others to propose strong evolution in galaxy luminosities, 
perhaps with the rate of evolution itself a function of luminosity (Broadhurst et al. 1988; 
Babul & Rees 1992; Babul & Ferguson 1996). Still another approach, taken by [Gronwall 
fc Koo (1995)[ , is to derive local luminosity functions by hnding functions that can explain 
as well as possible the faint galaxy number counts without invoking strong evolution. The 
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Figure 8: Cumulative luminosity densities for LSB galaxies and for galaxies from the 
CfA redshift survey, plotted linearly as a function of Log L. The top axis indicates the 
corresponding B magnitude. The points represent the binned model independent LF for 
the LSB galaxies, the solid represents the maximum likelihood Schechter function for the 
LSB galaxies, and the long-dashed line represents the two-component model LF for the LSB 
galaxies. The short-dashed line is the Schechter function for all morphological types in the 
CfA redshift survey determined by Marzke et al. (1994) and the dotted line is the Schechter 
function determined by Marzke et al. (1994) for irregulars. The cumulation runs from high 
to low luminosities (i.e., from right to left). 


luminosity functions they derive predict more local low-luminosity galaxies than are observed 
in existing surveys. At the very least, recent surveys for LSB galaxies indicate that the 
galaxy density at z = 0 is higher than previously assumed which means, at some level, the 
apparent excess of faint galaxies at high redshifts is at least in part an artifact of improper 
normalization at z = 0. The issue is how large this effect really is. 

McGaugh (1994) suggested that LSB galaxies such as those in the present sample could 
help reconcile the differences between observed local populations and this population of faint 
blue galaxies (hereafter, FBGs). He noted that, like the FBGs, LSB galaxies are generally 
blue (McGaugh & Bothun 1994) and weakly clustered (Mo et al. 1994). Furthermore, if 
current models of slow, continuous star formation LSB galaxies are correct (McGaugh & 
Bothun 1994), ^cGaugh (1994)| argued that LSB galaxies should become only slightly redder 
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over the timescales of interest, 0 < z ^ 0.5. He also demonstrated through a simple analytic 
calculation that the deep CCD surveys would be more sensitive to LSB galaxies at z ^ 0.4 
than wide-held photographic surveys are to local {z ^ 0.1) LSB galaxies. He argued that 
including nearby LSB galaxies in the local luminosity function could reconcile the number 
of low-luminosity galaxies in the local population with the FBG population. 


More recently still, priver et ah (1995b)| and [Driver et ah (1995a)| have examined 
the morphological mix of the faint held galaxies using data from the Hubble Space 
Telescope, down to a hux limit of m/ = 24.25. They compared the observed diherential 
number counts (number per square degree as a function of apparent magnitude) for three 
diherent morphological groupings to the predictions of various models. They concluded 
that diherential number counts of ellipticals and early-type spirals are consistent with the 
predictions of a no-evolution model based on standard local LFs (after renormalization) for 
these types taken from [Marzke et al. (1994a )| and [Loveday et ah (1992)| . They also found that 
the observed number counts of late-type spirals and irregulars were substantially in excess 
of similar no-evolution predictions for these classes. They could reconcile prediction with 
observation for this morphological class only by including a substantial amount of luminosity 
evolution (1.3 magnitudes of brightening by 2 ; ~ 0.5 for an Irr LF from Marzke et al. (1994)) 
or by using a sharply increased normalization (0* = 3.5 x 10 “^/i^qq Mpc“^; compare with 
Table |^) for the Irr LF. 


Our results underscore the uncertainty in the faint end slope of the held galaxy 
luminosity function. Priver fc Phillipps (199611 have shown that existing wide held redshift 
surveys place few constraints on the shape of the luminosity function below Mb = —16.5 
assuming Hq = 100 hioo kms“^ Mpc“^. By reaching lower in surface brightness, we have 
isolated a population of blue, LSB dwarfs that is absent from published luminosity functions, 
and which contributes strongly below Mb = —16. [Marzke et al. (1994a)| also saw evidence 
for a sharp upturn in the dwarf and irregular population at about the same luminosity. 
More recently, [Zucca fc et al. (1996)[ have shown that the luminosity function of [Lovedayj 
et al. (1992)[ calculated from the Stromlo-APM survey is signihcantly incomplete. The new 
ESO Slice luminosity function shows an upturn below Mb = —16, due to blue, star-forming 
galaxies, made up of a mixture of compact dwarfs and LSB galaxies. In all these studies, 
as in the deeper HST surveys of priver et al. (1995a)[ , the data are well described by a 
hybrid luminosity function consisting of a bright end Schechter function with a = — 1, and 
a faint end {Mb ^ —16 assuming Ho = 100 hioo kms“^ Mpc“^) power law with a slope 
—1.4 < a < —1.8. These faint galaxies are not major contributors to the luminosity density 
of the universe, but because the trend of M/L with luminosity is not well understood their 
contribution to the mass density is an open issue. 


There is also now evidence that some evolution may have occurred in the late-type 
galaxy population over the range 0 < ^ < 1. Lilly et al. (1995J have studied the evolution of 
the LFs over this range using date from the recent Canada-France Redshift Survey (CFRS). 
They found that the LF for red galaxies shows little change in number density or luminosity 
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over this range in z, but that the LF for blue galaxies appears to have brightened uniformly 
by about 1 magnitude by z ~ 0.75. 


Thus, the “excess” of FBGs consists of late-type spirals and irregulars, the same types 
which dominate the population of LSB galaxies found by the APM survey (Fig|^). These 
LSB galaxies have a Schechter function normalization approximately 6 times as large as that 
found by |Marzke et ah (1994a)| for HSB irregulars, so they expand the known local population 
well beyond that used by Priver et ah (1995b)| in their modeling. The LSB normalization 
still is not as large as that found necessary by priver et al. (1995b)| to account for the FBG 
counts with no evolution, but the heretofore-uncounted LSB galaxies do help considerably to 
close the gap between local population estimates and the FBG counts. Taking this increase 
together with the modest evolution observed in blue galaxy LFs by |Lilly et al. (1995]], it may 
now be possible to make an essentially complete reconciliation between local populations and 
the FBGs. However, any such reconciliation will also require a model for the evolution of 
the FBGs which accounts for the very blue colors of both the FBGs and the local LSBs (see 
McGaugh & Bothun 1994, McGaugh 1994, and McGaugh et al. 1995a). 

We can demonstrate the rough equivalence between the local LSB dwarf population and 
the MDS population of Driver et al. (1995a) , which we presume to be at typical redshifts 
0-3 ^ ^ ^ 0.6 based on [Lilly et al. (1995)| . The LSB dwarfs with Mb > —16 are at typical 


distances of 10 ^ d ^ 40 hboo Mpc. They have central surface brightness in the range 
22.0 < /iB(0) ^ 25 magarcsec"^, and effective angular radii of 6 ^ r^ff ^ 20 arcseconds. 
If they are related to the LSB dwarfs in clusters, we expect them to have B — V ~ 0.5 
(Impey, Bothun, & Malin 1988). The late-type and irregular (Sdm/Irr) MDS galaxies have 
effective radii with a median value of 0.4 arcseconds (Im et al. 1995), which would scale 
to 20 ^ Te// ^ 40 arcseconds for a local population. The central surface brightnesses of 
disk-dominated categories in the MDS sample are 20 ^ p-/(0) ^ 22 magarcsec“^per [Mutz 
et al. (1994)1 , which is equivalent to 22 ^ /^s(0) ^ 24 magarcsec”^ locally, assuming no 
evolution. The Sdm/Irr galaxies have observed colors of V — I ~ 1 per [Gasertano et al. 
( 1995)[ consistent with a local star-forming dwarf color oi B — V ~ 0.5, once again assuming 


no evolution. There could well be a subset of the MDS population which fade in surface 
brightness and redden to below the detection threshold for our blue photographic survey. 
Such galaxies would be absent from all local catalogs. 


6 Conclusions 

We have estimated a luminosity function for galaxies with surface brightnesses fainter than 
/i(0) = 22.0 magarcsec“^, which is the approximate faint limit of /i(0) for galaxies covered 
by the GfA Redshift Survey. We hnd that this LSB LF has a steeply rising tail at low 
luminosities (a = 1.42), comparable to that found by [Marzke et al. (1994a)| for galaxy types 
8 < T < 10. The LSB LF has a normalization lower than that found for the overall GfA 
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survey, but much higher than that found for types 8 < T < 10. Thus estimates of the total 
population of local galaxies based on the CfA survey are missing at least one-third of the total 
number of galaxies due to surface brightness selection bias. These previously unaccounted- 
for LSB galaxies can help considerably to resolve the apparent difference between estimates 
of the local population and the large numbers of faint blue galaxies observed at moderate 
redshift. 
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